A D6J8749 


ORC  65-15 
JUNV  1965 


im 


m 


UNIVERSITY  OF  CALIFORNIA  -  BERKELEY 


SEMI -NARROW  CHANNELS 


by 


James  Bigelow 
Operations  Research  Center 
University  of  California,  Berkeley 


June  1965 


ORC  65-15 


This  research  has  been  partially  supported  by  the  Office  of  Naval  Research 
under  contract  Nonr-3656(02 ),  and  with  the  United  Spates  Public  Health 
Service  under  contract  GM-96O6  with  the  University  of  California.  Repro¬ 
duction  in  whole  or  in  part  is  permitted  for  any  purpose  of  the  United 
Stales  Government. 


SEMI -NARROW  CHANNELS 
by 

James  Bigelow 
• 

Consider  a  channel  or  passage  connecting  compartments  A  and.  B  that 
is  wide  enough  so  that  two  particles  in  the  channel  can  pass  one  another. 

It  is  the  purpose  of  this  paper  to  derive  the  laws  giveming  the  flow  of  par¬ 
ticles  in  either  direction  through  such  a  channel. 

Side  A  *a  8  b  &  »  b  b  H  Side  B 
a  b 

Hie  position  of  a  particle  in  the  channel  can  be  found  by  projecting 
the  centers  of  all  the  particles  onto  the  axis  of  the  channel.  If  a  parti¬ 
cular  point  is  k  from  the  left,  then  the  particle  to  whose  center  the  point 
corresponds  is  in  the  position.  Hie  assumption  is  that  the  centers  of 

any  two  particles  will  project  onto  different  points  of  the  axis. 

An  event  occurs  if  a  particle  enters  the  channel  from  either  side  A  or 
B  .  It  may  either  strike  a  given  particle  in  the  channel  with  probability 
r  ,  or  bypass  it  with  probability  s  »  1  -  r  .  If  it  bypasses  the  particle, 
it  will  continue  moving,  and  either  strike  or  bypass  the  next  particla  a- 
long  its  path.  If,  cn  the  other  hand,  it  strikes  (say)  particle  k  ,  the 
moving  particle  comes  to  rest  in  position  k  ,  and  the  struck  particle  moves 
in  the  direction  of  the  original  motion.  Hius,  if  a  particle  enters  from 
A  ,  any  particle  it  strikes  will  move  toward  B  .  Also,  once  a  particle  is 
moving,  it  follows  the  same  laws  of  motion  as  the  entering  particle. 

We  have  assumed  here  that  the  probability  r  ,  that  a  moving  particle 
strikes  a  stationary  particle  in  position  k  ,  is  independent  of  k  .  With 
or  without  this  assumption,  however,  when  an  event  occurs,  exactly  one  par¬ 
ticle  will  move  to  position  k  and  one  particle  will  continue  beyond  po- 
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sition  k  for  every  k  .  Hence,  the  motion  traverses  every  position  in 
the  channel,  and  any  particle  may  be  struck. 

If  we  adopt  these  laws  of  motion,  we  may  conclude  that  if  a  particle 
enters  from  A  ,  sooner  or  later  a  particle,  though  not  necessarily  the 
same  one,  will  exit  into  side  B  .  Hence,  the  channel  always  holds  the  same 
number  of  particles  n  .  We  let  n  be  the  length  of  the  channel. 

We  will  assume  that  simultaneous  entries  from  both  sides  do  not  occur. 
Similarly,  we  shall  say  that  once  a  particle  enters  the  channel,  it  is  such 
a  short  time  before  some  particle  leaves  the  channel  as  a  result,  that  no 
other  entry  occurs  before  this  happens.  This  assumption  is  stronger  than 
necessary,  but  it  has  the  advantage  of  simplicity. 

Finally,  we  let  'p'  be  the  probability  that  the  next  entry  will  be 
from  side  A  ,  and  q  ■  1  -  p  be  the  probability  that  the  next  entry  will  be 
from  side  B  .  Then  we  can  calculate  the  probability  that  a  particle  in  po¬ 
sition  k  will  move  toward  (say)  side  B  .  as  a  result  of  an  entry.  Clearly 
this  requires  an  entry  from  side  A  ,  and  that  some  combination  of  particles 
hitting  and  bypassing  others  results  in  the  particle  ch  is  moving 
striking  particle  k  ,  when  the  motion  reaches  position  k  .  For  k  ■  1  , 
clearly, 


pr (particle  1  moves  toward  B)  ■  pr 

For  k  -  2  ,  the  entering  particle  may  miss  particle  1  and  hit  2,  with  prob  ■ 

? 

prs  ,  or  hit  1,  which  then  strikes  2,  with  prob  ■  pr  .  Thus, 

2 

pr  (particle  2  moves  toward  B)  ■  prs  +  pr 

■  pr  ( 6  *■  V ) 

n  pr 


Similarly,  for  k  ■  3  ,  the  entering  particle  can  miss  1  and  2,  and  hit 
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2 

3,  with  prob  ■  prs  ,  or  miss  1  and  hit  2,  which  then  strikes  3>  with  prob  ■ 

2  2 
prs,  or  it  can  hit  1,  which  then  may  miss  2  and  hit  3,  with  prob  ■  pr  s,  or 

the  entering  particle  may  hit  1,  which  may  hit  2,  which  in  turn  strikes  3, 

3 

with  prob  ■  prJ  .  Then: 

2  2^ 

pr(particle  3  moves  toward  B)  ■  prs  +  2pr  s  +  pr3 

■  pr(s^  +  2rs  +  r^ ) 

2 

-  pr  (s+r) 

-  pr 

k-1 

In  general,  there  are  (  j  )  ways  in  which  J  of  the  first  k-1  par¬ 
ticles  may  move  while  the  others  stand  still.  Hence,  the  probability  that 
this  occurs  is  Just: 

pr(exactly  J  of  the  first  k-1  are  struck)  ■  (  j^Jr^s^  ^  1 

and  the  probability  that  the  k^  will  be  struck  is  Just: 

k-1 

pr(k  will  move  toward  B)  ■  pr  ^  (^j^)r^sk  ^  ^ 

J-0 

■  pr(r+s  )k  ^ 

-  pr 

One  can  argue  this  another  way.  Since  r  is  the  probability  that  the 
moving  particle  will  hit  a  given  particle,  r  should  be  the  probability 
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that  the  given  particle  will  be  hit  by  the  moving  particle.  Thus  pr  is 
the  probability  that  it  ib  hit  from  A  ,  and  hence  moves  toward  B  .  Simi¬ 
larly,  the  probability  that  particle  k  will  move  toward  A  is  qr  .  Then 
the  probability  that  k  will  move  is  pr  +  qr  ■  r  ,  and  the  probability 
that  it  remains  stationary  is  a  result  of  an  entry  is  1-r  ■  s  . 

Now  it  is  easy  to  calculate  the  probability  that  particle  k  ,  once  it 
is  moving  toward  B  ,  will  Jump  to  position  J  ,  for  J  >  k  .  It  must  miss 
all  the  (j-k-l)  particles  between,  with  probability  s^"^-1  ,  and  strike 
J  .  T^ius, 

8J*k*1r  if  J  >  k 

pr(k  -*  J  |  k  moves  toward  B)  ■  j  b  ifj*k 

^  0  otherwise 

sk‘,3‘1r  if  J  <  k 

and  pr(k  -*  J  |  k  moves  toward  A)*"js  ifj*k 

^  0  otherwise 

Hence,  we  can  easily  calculate  the  unconditional  probability  that  k  will 

move  to  J  as  a  result  of  the  next  entry. 

pr(k  ”♦  j)  *  pr(k  -*  J  |  k  moves  bowar§  B)  <■  p»(lrm»ts ‘toward  B) 

+  p(k  -*  J  |  k  moves  toward  A)  •  pr(k  moves  toward  A) 

pr2sJ"k“1  if  J  >  k 

Or  p(k  -*  J)  -  <|  qr2sk”1' -1  if  J  <  k 

s  if  J  ■  k 

We  list  below  all  the  possiole  movements  k  can  undergo  as  a  result  of  an 
entry  and  their  probabilities,  as  well  as  the  probability  of  movements 
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which  originate  or  end  at  either  side  A  or  B  . 

1)  For  a  particle  in  the  channel  at  position  k  to  move  to  position 
J  ,  as  a  result  of  the  next  entry: 

2  j-k-1  ^  . 

pr  sJ  ,  if  j  >  k 

P(k  -  J)  -^qr2sk“J_1  ,  if  J  <  k 

s  ,  if  J  -  k 

[ 

2)  For  a  particle  in  the  channel  at  position  k  to  move  to  either 
side : 

p(k  *♦  side  B)  ■  prs11  k 
p(k  -*■  side  A)  *  qrsk  1 

3)  For  a  particle  entering  the  channel  to  come  to  rest  at  position 
k  : 

p(side  A  •*  k)  ■  prsk  1 
p(side  B  -*  k)  ■  qrs11  k 

U)  For  a  particle  to  move  from  one  side  to  the  other,  never  striking 
a  particle  in  the  channel : 

p(fide  A  -*  side  B)  ■  psD 
p(side  B  -*■  side  A)  ■  qs11 

What  we  wish  to  find  is  the  probability  f^  that  a  particle  will  move 

from  side  A  to  B  given  an  entry,  and  the  probability  f  of  a  movement 

B 
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from  B  to  A  . 


where  y^  ■  the  probability  of  finding  an  'a'  particle  in  position  k  , 

yk  ■  1-y^  ■  the  probability  of  finding  a  *b '  particle  in  position  k  . 

Thus,  we  must  first  find  an  expression  for  the  y's  and  y's  .  However, 
if  the  channel  is  in  steady  state,  this  is  not  hard  to  do. 


fo\  k-1  A  V  2  k’J“1  A  \  2  J-k-1 

(2)  yR  -  syk  +  prs  ♦  ^  pr  8  yj  +  qr  8  yj 

J-l  J-k+1 

Fortunately,  it  is  unnecessary  to  solve  for  all  of  the  y^'s  explicitly. 

n 

_  \  n-k  ,  n 

fA  ’  L  pre  yk +  ps 

k«l 

n-1 

V  n-k  .  .  n 

-  [_>  Prs  yk  +  pryn  +  ps 

k-1 

n-1 

_  .  .  n-1  .  V  2  n-k-1 

But  yn  -  syn  ♦prs  +  ^  pr  s  yk  . 

k-1 
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If  r)IO, 


yn  -  pB 


i-l  .  \  n-k-1 

♦  /  Prs  2 


k 


k-1 


n-1 


Or 


syn  -  ps  + 


V1  n-k 

1  PrS  J 


k-l 


So:  f A  -  (pr  +  s)yn  -  (l-qr)yQ 
(3) 

Similarly,  f  -  (qr  +  s)^  -  (1-pr)^ 


As  for  the  special  case  r  -  0  . 

Then  f  -  p  ,  f  -  q  (from  (l)  ),  and  the  equations  (2)  reduce  to: 
A  B 


Lemma  1. 


-  BtVffiJ -  ,  for  0  <  r  <[  1  , 


where  D  »  the  determinant 


Of 


1 

-qr 

n-2  f 

— qr  s  ••••*—  qr  s 

-pr 

1 

n-3 

-qr  . -qrs 

-prs 

-pr 

n-4 

1  -qrs 

t 

• 

• 

• 

• 

n-2 

n-3 

• 

• 

• 

n-4 

-prs 

-prs  J 

-prs  1 

T 

Note  that  thie  matrix,  post-multiplied  by  the  vector  (y^  ,  y2  ,  ...,y  )  , 
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should  yield  p(l,  s,  s^, . . . ,  s11  ^  )^  ,  if  equations  (2)  are  to  be  satisfied 


for  r  ^  0  . 

Proof:  Following  Cramer’s  rule, 

1 

-pr 

Vn  "  _prs 


-qr 

1 

-pr 


-qrs 

-qr 


•  •  • 


P 

ps 

r 

C 

ps 


n-2  *  n-3 

-prs  -prs 


ps 


n-1 


We  can  move  the  last  column  (which  is  Just  the  right-hand  side  of  the  prob¬ 
lem  noted)  to  the  first  position.  Since  it  will  pass  (n-l)  columns,  this 
will  change  the  sign  by  (-l)n“^  . 


(■Dn\yn 

■  n  n 


P 

ps 

ps‘ 


ps 


n-l 


1 

-pr 

-prs 

r 

-prs 


-qr  qrs 

1 


n-3 


-pr 


n-2 


-pr 

1 

-pr 


By  successively  adding  (-s)  times  a  row  to  the  row  below  it  (starting  with 
row  n-l)  ,  we  achieve  an  upper- triangular  structure. 


(-s-pr) 


Garbage 


(-s-pr) 


0 

(-s-pr) 

What  the  garbage  is  does  not  matter,  for  the  determinant  of  the  matrix  is 
Just  the  product  of  its  diagonal  terms. 


-8- 


Hence, 


(-l)n"XDnyn  -  p(-s-pr  )n”1  -  p(l-qr)n'1(-l)n"1 


,,  \  p(i-qr)n_1 

(M  yn  -  - 

n 

(The  restriction  on  r  is  necessary  to  put  equations  (2)  into  the  ap¬ 
propriate  fOTl.  ) 

Corollary 


(5) 


qU-pr.y 


n-1 


n 


This  follows  from  the  symmetry  of  the  problem. 


Lemma  2.  Consider  as  a  function  of  p  .  Then 

D  (p)  -  <i(l-Pr),n  for  ^  1/2 

n  p-q  T  ' 

D  (1/2)  -  lim  D  (p) 
n  P  -1/2  n 

Proof:  Dn(p)  is  a  determinant.  Hence,  it  is  defined  as  a  polynomial  in 

p  ,  in  this  case,  considering  r  fixed.  Thus,  as  defined,  DQ(p)  is  con¬ 
tinuous  in  p  . 

If  DQ(p)  as  derived  ■  Dn (p )  as  defined,  for  all  p  ^  l/2  ,  then 

lim  D„(p)  as  derived  is  the  same  as  lim  D  (p)  as  defined  =  D  (l/2) 
P-l/2n  p-l/2n 

as  defined. 
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Hence,  if  we  ©*tt  show  that: 


we  are  Justified  in  letting 

D  (1/2)  -  lira  D  (p) 
n  P-  1/2  n 

To  show  (6),  look  at  Dq+1 

,  n-1 

1  -qr  qrs  -qrs 

-pr  1  qr 

D  .  -  -prs  -pr  1 

n+l 

-qr 

-prs11”1  -pr  1 

Adding  (-s)x  the  second  row  to  the  first  and  expanding  on  the  new  first 
row,  we  get: 

Dn+1  -  (1  +  prs)  Dn  ♦  (s  +  qr)  (-r)  (-l)n 

where  y|n^  ■  the  probability  of  finding  an  'a'  in  position  one  of  a  tun¬ 
nel  of  length  n  . 

Letting  y|n^"  1  -  y-[n^  >  anci  substituting  from  the  corollary  of  Lem¬ 
ma  1, 

Dn+1  ■  (l-qr )  Dn  +  qr(l-pr)n 

We  could  have  added  (-s)x  the  second  column  to  the  first  and  expanded  on 
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the  new  first  column  just  as  easily.  This  would  yield: 


Dn+1  »  (1-pr )  Dn  ♦  pr(l-qr)n 

For  p  /  l/2  ,  these  two  expressions  can  be  solved  for  . 

D  ,  pr(l-gr)n  -  pr(l-pr)11 

n  (1-qr)  -  (1-pr) 

or 

D  „  p(l-qr)n  -  q(l-pr)n 
n  p-q 

(This  proof  is  due  to  Dr.  S.  Levin.  An  alternate  inductive  proof  is  also 
possible.)  Q.E.D. 

Now  e  can  prove: 

Theorem  1  (f»  -  f„)  *  p-q  ,  independent  of  n  and  r  . 

A  D 

Proof:  Trivial  for  r  ■  0  .  since  f .  ■  p  ,  and  f_  «  q  . 

-  A  B 

(3)  For  r  /  0  ,  f A  *  (l-qr)yn 

fB  ■  (l-pr)^ 


f  -  f  .  (l-qr)  .  (l-pr)q(I-pr)n~1 

A  B  K  *  '  D  D 

n  n 

=  p(l-qr)n  -  q(l-pr)n 

,p(l-qr)n  -  q(l-pr)nx 
V  p-q 
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Q.E.D. 


Perhaps  more  interesting  is 


Theorem  2: 


p(i-qr )n 
q(l-pr  )n 


for  0  £  r  £  1  . 


Proof :  At  r 


from  the  special  case  mentioned  earlier.  At 


any  other  r  , 


(8) 


(l-qr)yn  p(l-qr)n/Dn 
h  '  =  q(l-pr)n/Dn 

=  p(l-qr)n 

fB  q(l-pr )n  Q.E.D. 


f^  and  f  are  merely  probabilities  that  flow  will  occur  in  one  di¬ 
rection  or  the  other  as  a  result  of  a  single  entry.  To  calculate  the  ac¬ 
tual  rate  of  flow  we  must  replace  the  probability  of  an  entry  from  either 
side  with  the  rate  at  which  entries  from  either  side  occur.  For  a  chemical 
solution,  we  allow  these  rates  on  either  side  to  be  proportional  to  the 
concentration  of  the  entering  particles  on  that  side,  but  not  necessarily 
with  the  same  proportionality  factor.  Differences  could  occur  because  of  a 
temperature  difference,  a  pressure  difference,  or  an  electric  potential  dif¬ 
ference  between  the  ends  of  the  channel. 

Also,  since  many  different  types  of  particles  might  enter  from  either 
side,  we  must  let  the  entry  rate  from  each  side  to  be  a  weighted  sum  of  the 
concentrations  of  the  entering  particles.  The  weights  will  be  called  entry 
constants,  and  designated  by  ^  (for  side  A)  and  ^  (for  side  B)  . 
Similarly,  »  concentration  of  the  ith  entering  species  on  side  A  , 
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*  concentration  on  side  B 


Then 


P  "  k  r  Mi°i 


q  =  k 


I  /*»  I 

,  vici 


Now  define 


2  u,  C 


i  i 


2  u'Ci 


8 


1-qr  X  +  P 
1-pr  =  1  +  Xs 


Then : 


(9) 


fa  *  (z  ^ici 


n  1-X 


i-8nx 


(2  n’Cp 


1-X 

i-8nx 


Let  us  examine  the  theory. 


I.  Limiting  Cases 

As  s  -»  0  we  should  have  a  tunnel.  But,  lim  8 

s  -*•  0 


FA  - J  ^ici  ^rrsw’ 

A, 


h‘z^cl  (77^i) 

-A.  "*  A. 


X  .  Hence, 
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So: 


(10) 


n  Z  ^iCi 


n+1 


and  Fa  -  Fb  -  (Z  UiC1  -  Z  n^C’) 

These  two  equations  characterize  a  tunnel.  Hence,  this  limiting  case  is  va 
lid. 

As  s  -*■  1  we  should  approximate  simple  diffusion.  Note  that 

lim  3*1. 
s  *♦  1 

Thus,  FA  "  2  MiCi 

fb  ■ 

Hence, 

F. 

(11) 

fb 


F.  -  F 


B 


MiCi 


-  Z 


»*iCJ 


Again,  these  equations  characterize  simple  diffusion. 


II.  Approximation 

We  c£n  try  to  approximate  a  semi -narrow  channel  by  a  composite  channel; 

that  is,  a  weighted  average  of  a  free  diffusion  path  and  a  narrow  tunnel. 

We  can  assign  weighted  t)  to  the  tunnel  and  (1-tj)  to  the  diffusion  path. 
TP  S 

If  F  >  F  ,  and  denote  flow  through  the  tunnel,  diffusion  path,  and  semi¬ 
narrow  channel,  we  wish: 
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n  *  (1-1)  rj  - 

iFB*  FB 

Substituting: 


where  the  n  *  the  length  of  the  semi -narrow  channel,  which  is  known,  and 
m  ■  the  length  of  the  tunnel,  which  is  to  be  determined. 


For  an  exact  solution,  we  must  have  ■  T)  .  This  condition  is  satis¬ 
fied  if  we  let: 

D  S 
U  ■  u 

„  .  T  S 

and  ■  n 


Then  T)  and  t]  reduce  to: 


1  -  t\ 


-15- 


Think  of  T]  &b  a  function  of  X  ,  and  let  X  0  . 


Then 


T)(0 )  -  1  -  Sn  ,  since  3  -  rTTs 
If  tj  is  to  be  constant,  then: 

(i-en)0--xm+1)  -  (l-dn)(l-pnX)(l^“) 

Or,  substituting  for  3  , 

t(l+XS)D  -  (X+sfl  t'.-x"*1]  -  (i-sn)t(i+xs)n  -  x(x  s)n](i-xn)  . 


Comparing  constant  and  first  order  coefficients  of  both  polynomials,  we  get 
that  these  relations  must  hold: 


(1-Sn)  ■  (1-Sn)  (constant  terms) 


and 


(nXS  -  nSn”1x)  «  (l-Sn)fnXS  -  XSn]  (first-order  terms) 


Or 


n-2  n-1  x  n  „2 n-1 

S  -  S  +  S  -  S 


This  is  clearly  untrue.  Thus,  a  semi -narrow  channel  is  not  a  composite. 
However,  if  we  examine  r\  more,  we  see  that 


i  -  en  -  xm+1  ♦ 

|0  —■■■*■ . . 

1  «  m  ,  «n.  m+1 


m+1 


l  -  p  x  -  x“  +  3  x 


Restricting  t]  £  1 


i  -  en  -  x"**1  +  pnx"M'1  ^  i  -  A  -  xm  ♦  p“x“  1 
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Or 

$n(\  -  1)  £  -  1) 

Thus, 

3n  £  x“  for  X  >  1 

0n  ;>  X®  for  X  <  1 

If  we  graph  the  two  functions  of  X  ,  that  is,  pn  and  xm  ,  we  see 
the  following: 


then  we  have  T)  ■  1  ,  and  that  for  X  near  1  ,  the  semi -narrow  channel  of 
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length  n  acts  like  a  tunnel  of  length  m  .  We  find  the  value  m  by  equa¬ 
ting  the  derivatives  at  X  *  1  .  This  is  a  first-order  approximation. 


Now 


np 


n-1  dfl  ,m-l 
-  mX 


dX 


P- 


X+S 

1+XS 


d £m  1-S2 
dx  (1+XS  )2 


Evaluating  at  X  ■  1  , 


m  *  n 

Notice  that  as  n  gets  very  large,  the  approximation  is  valid  over  a 
larger  and  larger  range.  In  fact,  for  n  -*■  »  ,  the  semi-narrow  channel  be¬ 
comes  a  tunnel,  Just  as  for  n  -*■  0  ,  the  semi -narrow  channel  is  the  same  as 
a  diffusion  path. 


III.  Flows  Against  the  Gradient 


Notice  that  the  flow  of  the  i^*1  entering  species  is  Just: 


•i  -  {  -i  °i  en  -  }  ! 


1  -  X 

1  -  p“x  J 


Even  if  n  provided  p  >  1  ,  ard  n  is  large  enough,  i.e.,  if 


^ 


C1  "1  a 

,  or  if  —  <  — 7  p  we  can  have  species  i  flow  against  its 
uiCi  Ci  ui 


gradient.  Note  here  that  when  X  >  1  ,  then  P  >  1  ,  and  P  <  X,  so  that 
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a  tunnel  of  length  n  will  transport  species  i  agt.inst  a  higher  gradient 
than  a  semi -narrow  channel. 

IV.  Diffusion  vs.  Osmotic  Diffusion  Coefficients 

In  the  hunan  red  blood  cell,  measurements  of  permeability  to  water  are 
made  in  two  ways.  If  an  actual  osmotic  pressure  is  introduced,  the  net  rate 
of  flow  of  H^O  can  be  measured  and  an  osmotic  coefficient  calculated.  If 
osmotic  balance  is  achieved,  a  diffusion  coefficient  can  be  calculated  from 
the  unidirectional  water  flow  as  measured  with  labeled  water,  or  DHO.  The 
surprising  fact  is  that: 


4  “  2.5  4 , . 

osm  dlff. 

Such  a  difference  is  readily  explained  by  this  theory.  Let  us  propose 
channels  of  a  semi-narrow  type  through  which  only  water  flows,  to  a  first  and 
probably  good  approximation.  Then  we  assume  equal  temperature  and  pressure 
on  both  sides  of  these  channels,  and  we  say  that  water  is  unaffected  by  the 
membrane  potential.  Thus,  the  entry  constants  are  the  same  from  each  side. 
Finally,  we  treat  DHO  as  if  the  channel  cannot  distinguish  it  from  H^O  . 

If  there  is  an  osmotic  gradient: 

fa  -  fb  ■  u(c-c') 

Thus, 

4  “4 

osm 

Now,  if  osmotic  equilibrium  is  the  case,  then  X.  -  p  ■  1  .  But: 
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Using  L'hopital's  Rule:  *  - - - 


Clearly,  for  s  <  1  ,  ^  <  Mosm  •  If  we  put  in  n  fibers: 


To  calculate  r  ,  let: 


d  ■  the  diameter  of  an  entering  particle,  where  we 
assume  all  particles  are  of  equal  size. 

D  ■  diameter  of  the  channel,  letting  both  be  circular. 
Then  to  enter  the  channel,  the  center  of  the  particle  must  be  somewhere  in 


If  the  particle  is  to  strike  another,  the  centers  of  the  two  particles  must 
be  within  d  of  one  another,  giving  a  "striking  area"  of 


nd 

"  ~T~ 
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,  when  this  is  less  than  one.  Hence, 


We  let  r 


d  ,  ,  D 

- 2  »  for  d  <  2 


(D-d) 


1  ,  for  D  ^  § 


Now,  the  size  of  a  water  molecule  is  said  to  be  about  3&  ,  and  the 
length  of  the  channel,  which  we  assume  is  Just  the  thickness  of  the  membrane, 
is  10QA  .  Assuming  close  packing  of  particles  in  the  channel  (i.e.,  each 
molecule  requires  only  3&  of  length). 


n 


100 

3 


33 


Thus,  if 


s  -  .914 

or  r  ■  .086 

0 

This  gives  D  ■  13.2A 

This  model  thus  predicts  that  water  channels  will  have  a  diameter  of 

O 

13. 2A  .  Solomon  has  indirectly  measured  the  channels  and  calculates  that 

0 

their  diameter  is  about  10A  .  We  are,  then,  at  least  in  the  right  ball  park. 
If  we  do  not  assume  close  packing,  and  instead  let  fewer  water  molecules  re¬ 
side  in  the  channel,  then  the  channel  diameter  is  smaller.  Thus,  I3.2A  is 

an  upper  bound.  Similarly, - may  be  greater  than  2.5  .  This  is  the 

^diff 
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smallest  measurement  in  the  literature.  As  the  ratio  grows,  the  channel  will 


shrink. 

Although  the  mathematics  are  ironclad,  this  model  of  flow  through  semi- 
narrow  channels  is  only  an  approximation  to  what  may  really  happen.  Some  of 
its  weaknesses  are  obvious.  For  example,  we  assumed  that  the  probability  r 
of  a  hit  was  independent  of  the  position  in  the  channel  which  the  struck  par¬ 
ticle  occupied.  One  might  relax  this  assumption  and  subscribe  to  the  belief 
that  parti cles  near  the  end  of  the  channel  should  be  easier  (or  harder )  to 
hit  than  those  in  the  center.  One  reason  that  particles  near  the  end  of  the 
channel  might  escape  being  struck  as  often  as  those  in  the  center  could  be 
that  center  particles  could  be  packed  more  densely  than  the  end  particles, 
and  hence  fill  up  a  greater  part  of  the  cross-sectional  area  of  the  center 
of  the  channel.  In  the  figure,  one  can  6ee  how  this  might  operate. 

O  "" 

O— *  o  o 

The  moving  particle  would  be  certain  to  strike  one  of  the  particles. 

A  second  weakness  of  the  model  is  that  it  assumes  that  any  particle 
which  enters  the  channel  has  the  same  chance  of  hitting  a  given  particle 
k  as  any  other  enterings  particle  ,  whereas  one  entering  particle  may  be  of 
a  different  size  than  another.  If,  for  example,  we  suppose  that  two  parti¬ 
cles  of  different  sizes  are  permitted  in  the  channel  with  diameters  d^  <  , 

a  particle  of  type  1  would  be  less  likely  to  strike  another  of  type  1  than  of 
type  2  .  Similarly,  a  particle  of  type  2  would  be  less  likely  to  strike  one 
of  the  first  type  than  the  second.  In  fact,  a  channel  could  act  like  a  tunnel 
to  the  larger  particle,  and  a  semi-narrow  channel  to  the  smaller.  If,  for 

example,  d?  -  ^  +  c  ,  and  d^  «  ^  ,  then  we  might  have  a  channel  operating  as 
a  tunnel  for  large  particles  (type  2)  and  as  a  free  diffusion  channel  for  the 
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small  (type  1)  particles.  Alternately,  a  large  particle  might  act  as  an 
effective  plug  to  smaller  ones,  and  hence  slow  down  or  stop  what  would  other¬ 
wise  be  nearly  free  diffusion.  In  short,  this  model  can  and  should  be  re¬ 
fined. 
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